New Development of Monte Carlo Techniques for Studying Bottle-brush Polymers 
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Due to the complex characteristics of bottle-brush polymers, it became a challenge to develop 
an efficient algorithm for studying such macromolecules under various solvent conditions or some 
constraints in the space by using computer simulations. In the limit of a bottle-brush polymer 
with a rather stiff backbone (straight rigid backbone) , we generalize the variant of the biased chain 
growth algorithm, the pruned-enriched Rosenbluth method, for simulating polymers with complex 
architecture, from star polymers to bottle-brush polymers, on the simple cubic lattice. With the 
high statistics of our Monte Carlo results, we check the theoretical predictions of side chain behavior 
and radial monomer density profile. For the comparison of the experimental data for bottle-brush 
polymers with a flexible backbone and flexible side chains, based on the bond fluctuation model 
we propose another fast Monte Carlo algorithm combining the local moves, the pivot move, and an 
adjustable simulation lattice box. By monitoring the autocorrelation functions of gyration radii for 
the side chains and for the backbone, we see that for fixed side chain length there is no change in 
the behavior of these two functions as the backbone length increases. Our extensive results cover 
the range which is accessible for the comparison to experimental data and for the checking of the 
theoretically predicted scaling laws. 



I. INTRODUCTION 

The so-called "Bottle-Brush" polymers consist of a 
long molecule serving as a "backbone" on which many 
side chains are densely grafted. The conformational 
change of bottle-brush polymers is mainly caused by the 
following factors: backbone length, side chain length, 
grafting density, type of monomers (chemical compound), 
and solvent quality which can be adjusted by changing 
the temperature, pH value, etc. In the previous Monte 
Carlo studies of bottle-brush polymers, both coarse- 
grained models on lattice [T]-[3] and on off-lattice [4]-[6], 
show that it is difficult to obtain high accuracy results for 
simulating large bottle-brush polymers with high graft- 
ing densities. 

On a coarse-grained scale, the bottle-brush polymer 
with densely grafted side chains may resemble a flexible 
long sphero-cylinder [T, 8 . The complicated structure of 
bottle-brush polymers is therefore described in terms of 
multi-length scales such as the contour length Lcc, the 
end-to-end distance of the backbone, Reh^ and of the side 
chain, R^^ the cross sectional radius Rcs^ and also the per- 
sistence length Ip which describes the intrinsic stiffness 
of the backbone, i.e., within the distance ^p, the cylinder 
is approximately straight. With computer simulations, 
one can estimate not only all these length scales but also 
those physical quantities measured by experiments such 
as the structure factors S{q) which describe the scatter- 
ing function from any part of the bottle-brush polymers, 
and the radial monomer density profile p(r). Therefore, 
it is necessary to develop an efficient algorithm for a 
deeper understanding of the complicated structures of 
bottle-brush polymers with larger size and higher graft- 
ing densities in order to control their functions for the 
applications in industry. 

In this article, we first explain the models and the 
algorithms for simulating bottle-brush polymers with a 
straight rigid backbone, and with a flexible backbone in 



Sec. [Hj For simplicity, here we only focus on the case 
that the bottle-brush polymers are under good solvent 
conditions. In Sec. |lll| we present our results and ex- 
plain the connections between these estimates obtained 
by computer simulations, the theoretical predictions, and 
the experimental data, respectively. Finally we give some 
conclusions in Sec. HVl 



II. MODELS AND ALGORITHMS 

For studying bottle-brush polymers under a good sol- 
vent condition, we first consider that bottle-brush poly- 
mers consist of a rigid backbone and flexible side chains. 
A simple coarse-grained model on a simple cubic lat- 
tice is used, where the backbone is simply a rigid rod 
and flexible side chains are described by self- avoiding 
walks (SAWs) so that no multi-occupation of monomers 
on the same site is allowed. We apply a biased chain 
growth algorithm with resampling which is a variant of 
the pruned-enriched Rosenbluth method (PERM) [9HTT] 
for the simulations. For the comparison between experi- 
mental data and Monte Carlo simulations of bottle-brush 
polymers, we need to consider a more complicated case 
that the backbone is also flexible. The bond fluctuation 
model [12-15 is used, where the backbone and all side 
chains are described by SAWs on a sim ple cubic lattice 
but with some constraints (see sec.|IIB). We propose an 



algorithm which combines the local 26 moves, the pivot 
moves and an adjustable simulation lattice box (LPB) 
for the simulations [71 H] [16] 



A. Simple coarse-grained lattice model with PERM 

In the simple coarse-grained model, the backbone is 
fixed on the simple cubic lattice in the direction along 
the z-axis. A/'^ monomers of the backbone are located on 
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FIG. 1: (a) Schematic drawing of the geometric arrangement for the simple coarse-grained model, (b) A snapshot of a bottle- 
brush polymer with Nb = 128, N = 2000, and a = 1/4 on the simple cubic lattice. Note that different colors are used in order 
to distinguish between different side chains, and the periodic boundary condition is undone for the sake of better visualization. 



the lattice sites as shown in Fig. [T] ric side chains con- 
sisting of N monomers each are grafted to the backbone 
monomers with equal distance 1/a between two succes- 
sive grafting sites on the backbone, where a is the graft- 
ing density defined by a = ric/Nij. Since mainly we want 
to check the scaling laws for very long side chains, the 
periodic boundary condition is introduced in the direc- 
tion along the backbone to avoid the end-effects associ- 
ated with a finite backbone length. Differently from the 
conventional MC method where an initial configuration 
is set up as a starter for the simulation, the conforma- 
tion of a bottle-brush polymer is built by growing all side 
chains simultaneously with PERM. The partition sum of 
bottle-brush polymers of Uc side chains of length A^ each. 



Znu, = 7,17 
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walks 



is therefore the total number of all possible configurations 
of Tie interactive self-avoiding random walks of steps N . 
Only the excluded volume effect is considered here. 

It is straight forward to apply the similar method for 
growing a star polymer [10] to the simulations for grow- 
ing a bottle-brush polymer [11]. Only in the latter case, 
side chains can either be attached to the same site or the 
different sites on the backbone depending on the graft- 
ing densities. In the process of growing a bottle-brush 
polymer, one has to be aware that both the interactions 
between monomers in the same side chain, and the inter- 
actions between monomers on different side chains have 
to be taken into account. If one side chain is grown en- 
tirely before the next side chain is started, it will lead 
to a completely "wrong" direction of generating the con- 
figurations of a bottle-brush. Therefore, one has to use 
the strategy that all side chains are grown simultane- 
ously. Namely, a monomer is added to each side chain 
step by step until all side chains having the same length, 
then that the next round of monomers is added. Af- 



ter we labelled all monomers by numbers, it goes back 
to the problem of growing a linear chains from the first 
monomer to the {ncN)th monomer [9 . Using PERM, 
the configurations of bottle-brush polymers are built by 
adding one monomer at each step, and each configura- 
tion carries its own weight. A wide range of probability 
distributions can be used for selecting one of the near- 
est neighbor free sites of each side chain end at the next 
step, but the efficiency of the algorithm depends on the 
choice of the distribution. For the current problem, the 
bias of growing side chains is used by giving higher prob- 
abilities in the direction where there are more free next 
neighbor sites and in the outward directions perpendic- 
ular to the backbone, where the second part of bias de- 
creases with the length of side chains and increases with 
the grafting density. The total weight Wm (jn = nuc) 
for a bottle-brush polymer of all side chains having the 
length n with an unbiased sampling is determined recur- 
sively by Wm = nr=i^fc = Wm-iWm- As the weight 
Wm is gained at the mth-step with a probability p^, one 
has to use Wm/Pm instead of Wm- By taking the average 
of all possible configurations, the partition sum 
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can be estimated directly, where Mm is the total number 
of configurations {a}. This is the main advantage of 
using PERM. For any observable Am-, the mean value is 
therefore, 
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In order to suppress the huge fluctuation of the prob- 
ability distribution and enrich those configurations with 
high weight the population control is made in the way of 
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FIG. 2: (a) Schematic drawing of applying local 6 (L6) moves and local 26 (L26) moves to a monomer on the site of a simple 
cubic lattice. In the "L6" moves, a monomer is tried to move to the nearest neighbor sites in the six directions, while in the 
"L26" moves, a monomer is not only tried to move to the nearest neighbor sites but also to the next nearest neighbor sites 
and the sites at the 8 corners which are in v^ lattice spacings away from the chosen monomer, (b) Two types of pivot moves 
applied to a randomly chosen monomer on the backbone, and to a randomly chosen monomer on a randomly chosen side chain. 




FIG. 3: The two observables, the radius of gyration of back- 
bone monomers, Rgb, and the space occupation of bottle- 
brush polymers, (A^c, A^, A^^), are indicated in the schematic 
drawing by taking a snapshot of a bottle-brush polymer with 
Nb = 131, Ar = 6, and cr = 1. 



pruning low weight configurations and cloning those con- 
figurations with high weight. Two thresholds W^ and 
W^ are introduced, 



^+ = G+^. 
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where Zm is the current estimate of the partition sum 
{Eq. Q}, and C+ and C- are constants of order unity. 
The optimal ratio between C+ and C- is found to be 
Cj^/C- ~ 10 in general. For our simulations, we use 
W^ = oo and W^ = for the first configuration hitting 
all side chains of length n. For the following configura- 
tions, we use W+ = CZ^(c^/co) and W~ = 0.15W+, 
here C = 3.0, and Cm is the total number of configura- 
tions of all side chains having length n. If the current 
weight Wm{(y) > W^ for the configuration a, one pro- 
duces two identical copies of this configuration, replaces 
their weight Wm{(^) by Wmio^)/"^- If Wm{(^) < W~, one 
calls a random number r where r G [0, 1]. If r > 1/2, 
the configuration is kept but the weight is replaced by 



2VFm(a), while the configuration is killed if r < 1/2. Oth- 
erwise, the configuration is kept with the weight Wm{o^)- 
A typical configuration of bottle-brush polymers un- 
der a good solvent condition generated by PERM is 
shown in Fig. [ijb). It consists of A^5 = 128 backbone 
monomers, N = 2000 side chain monomers in each side 
chain, and the grafting density is 1/4. The total number 
of monomers is TVtot = A^6 + NaN^ = 64128. So far, it 
is the largest bottle-brush polymers in the equilibrium 
state, generated by MC simulations [11^. 



B. Bond fluctuation model with LPB 

For studying bottle-brush polymers with a flexible 
backbone and flexible side chains, we generalize the bond 
fluctuation model for a linear polymer chain to that for a 
bottle-brush polymer. In the standard bond fluctuation 
model [TJHIS], a flexible polymer chain is described by 
a SAW on a simple cubic lattice with bond constraints. 
Each effective monomer blocks all 8 corners of an elemen- 
tary cube of the lattice from further occupation. Two 
successive monomers along a chain are connected by a 
bond vector chosen from the set {(±2,0,0), (±2, ±1,0), 
(±2,±1,±1), (±2,±2,±1), (±3,0,0), (±3, ±1,0)}, in- 
cluding also all permutations. The geometry of a bottle- 
brush polymer with N^ backbone monomers, and with ric 
side chains of length N is arranged in the way that side 
chains are added to the backbone chain at regular spacing 
1/(7 = Niy/ric^ and two additional monomers are added to 
each chain end of the backbone. Thus, the total number 
of monomers is A/tot = [(^c — l)/<^ ± 1]±2. For our simu- 
lations, one of the simplest ways to set up the initial con- 
figuration is to assume that the backbone and side chains 
all have rod-like structures. Placing the backbone along 
the 2:-direction and fixing the bond length between two 
successive backbone monomers to be 3, and randomly 
choosing the bond vector of each side chain from one of 
the allowed bond vectors including all permutation in the 
xy-plane but keeping the bond vectors within each side 
chain fixed, the required condition of bond constraints 



is satisfied and no further check is needed. In our algo- 
rithm, instead of trying to move a chosen monomer to the 
nearest neighbor sites named by the local 6 ( "L6" ) moves 
for the standard bond fluctuation model, we use the local 
26 ("L26") moves [17 where it is tried to move to the 26 
neighbor sites as shown in Fig. [2|a). The local move is 
only accepted if the selected site is empty and the bond 
length constraints are satisfied. In addition, two types 
of pivot moves are attempted. One is that a monomer 
is chosen randomly on the backbone and the short part 
of the bottle-brush polymer is transformed by randomly 
applying one of the 48 symmetry operations (no change; 
rotations by 90^ and 180^; reflections and inversions). 
The other is that a monomer is chosen randomly from 
all side chain monomers, and the part of side chain from 
the selected monomer to the free end of the side chain is 
transformed by one of the 48 symmetry operations. 

We first make some test runs for a small bottle-brush 
polymers with Ni, = 32, iV = 6, 12, 24, and 48, and 
cr = 1 in order to compare the efficiency between the 
"L26" moves and the "L26" moves + pivot moves. For 
any observable A, the performance of the algorithm is 
determined by the autocorrelation function c(A,t), 
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Results of c{A,t) for the mean square gyration radius 
of the backbone, A = i?^^, and of the side chains, A = 
i?^^ (taking the average of all side chains at the same 
MC step t) plotted against the number of MC steps t 
show that the "L26" + pivot algorithm is two orders 
of magnitude faster than the "L26" algorithm for fixed 
side chain length N (figures are not shown here). In the 
"L26" algorithm, one MC step consists of A^tot "L26" 
moves, i.e., each monomer is selected once for the local 
move. In the "L26" + pivot algorithm, one MC step 
consists of A^tot "L26" moves, k^ times pivot moves of 
the backbone and kc times pivot moves of side chains, k^ 
is chosen such that the acceptance ratio is about 40% or 
even larger, while kc is nc/4. 

Two observables chosen for monitoring the equilibrat- 
ing process as shown in Fig. [3] are the radius of gyration 
of the backbone monomers, Rgb{t)^ and the space oc- 
cupation of the bottle-brush polymer in the Cartesian 
coordinates, {Ax{t)^ Ay{t), Az{t))^ since the average con- 
formations of bottle-brush polymers in equilibrium must 
be isotropic. For simulating small bottle-brush polymers, 
we can simply set the three orthogonal length scales 
in the Cartesian coordinates having equal length, e.g. 
Lx = Ly = Lz = SNi) in the equilibrating process. As the 
size of bottle-brush polymer increases, we will meet the 
problem of setting the simulation lattice box in our sim- 
ulations due to the limitation of the computer memory. 
The maximum volume of the box is V = L^LyLz = 2^^ 
for those computers we can access. The solution for it 
is to adjust the simulation lattice box during the equi- 
librating process and separate the process into several 
stages. Let's take a bottle-brush polymer with N^ = 515, 
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FIG. 4: Time series of the reseated square gyration radii for 
the backbone monomers, R^ij(t)/N^^ (a) and of the spacing 
occupation of the whole bottle-brush polymers in the Carte- 
sian coordinates, {Ax{t), Ay{t), Az{t)) (b). Results shown 
here are for bottle-brush polymers with Nb — 515 backbone 
monomers, N — 12 side chain monomers, and the grafting 
density a — 1. 



A/" = 12, and <j = 1 as an example. The equilibrating pro- 
cess is separated into four stages as follows. 



stage 1: 1 < 7V^^ < 128, L^ = 1545, L^^ = L 
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Stage 2: 1 < 7V^^ < 256, L^^ = 1201, L^^ = L 
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tf^ = 262144 MC steps 

stage 3: 1 < 7V^^ < 513, L^^^ 
tf^ = 262144 MC steps 



851, Ly 



(3) 



Lx 



(3) 



(4) 



Stage 4: 1 < A^^^ < 513, L^^ 
1310720 MC steps 
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Here A^^ is the pivot point selected from the backbone 
monomers. One has to be aware that the pivot points 
which can be selected for applying the pivot moves are 
also limited due to the current set up of the simulation 
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FIG. 5: Snapshots of bottle-brush polymers with Nt = 515, N = 12, and a = 1 in the equilibrating process including four 
stages. 



lattice box. However, once the system is in equilibrium 
one has to allow all possible moves, i.e. (1 < N^ < 
(7V5 — 2)). At every new stage /c, the initial configu- 
ration is taken from the last stage, and the size of the 
simulation lattice box is decided by the final space oc- 
cupation of the bottle-brush polymers at the last stage. 



I.e. 
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and Lf^Ll'^L):^ < 2^\ Time 



LT^ > A 
with t = ^ 

series of Rgb{t) and {Ax{t)^ Ay{t), Az{t)) are shown in 
Fig. [4] The four stages are separated by the vertical 
green curves. Taking some snapshots of the conforma- 
tions of the bottle-brush polymers at the Monte Carlo 
steps indicated by the arrows in Fig. |4jb) one can see 
how the conformations of bottle-brush polymers change 
during the equilibrating process as shown in Fig. [5] At 
the beginning, backbone and side chains are in rod-like 
structures. At the end of the first stage, only a small part 
of the backbone is flexible. As more backbone monomers 
are relaxed, we see that the backbone and the side chains 
become more and more flexible step by step. Finally an 
equilibrium state is reached. It takes about 1.25 hours 
CPU time on an Intel 2.8 GHZ PC for such a bottle- 
brush polymer to reach the equilibrium state by choosing 
h = 40 and kc = 128. 

It is more time consuming when pivot moves are ap- 
plied to the simulations for larger bottle-brush polymers. 
In order to know how much the efficiency is slowing down 
as the backbone length increases, we compare the au- 
tocorrelation functions c(i?^^,t) and c(i?^^,t) for bottle- 
brush polymers with N^j = 35 and N^j = 515 for three 
different side chain lengths A/" = 6, 12, and 24. For both 



cases in Fig. |6ja), we see that the decay of the autocor- 
relation function for the side chain structure occurs on 
the same time scale. It is also true for the backbone as 
shown in Fig. p[b) that the autocorrelation functions are 
plotted against the number of pivot moves t/c^. Clearly, 
the structural relaxation time is longer as the side chain 
length N increases [I83 . 



III. RESULTS 

According to the cylindrical geometry of bottle-brush 
polymers with a rigid backbone, we extend the Daoud- 
Cotton blob picture for a star polymer to that for a 
bottle-brush polymer. The space is partitioned into blobs 
of non-uniform size and shape. The blobs are not spheres 
but rather ellipsoids. Based on this theory, the scaling 
law for side chains in the radial direction (height of the 
bottle-brush) is given by. 
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where u is the Flory exponent for 3D SAWs {u ^ 0.588). 
With the first part of simulations by using PERM, results 
of the mean square height of bottle-brush polymers for 
three choices of backbone length TV^ = 32, 64, and 128, 
and several choices of the grafting density a from 1/128 
to 1 are shown in Fig. [TFa). We see that those curves 
of the same grafting density a coincide with each other. 
Increasing the grafting density a enhances the stretching 
of side chains. Considering that in the mushroom regime 
(cr ^ 0), the height of bottle-brush polymers should be- 
havior as 3D SAWs, i.e. Rh{N,a ^ 0) -- A/"^, one can 
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FIG. 6: Autocorrelation functions of the mean square gy- 
ration radii for the side chains c{R^cit) (taking the average 
of all side chains at each t) plotted against the number of 
Monte Carlo steps t (a), and for the backbone c{R'^i^, t) plotted 
against the number of pivot moves applied to the backbone, 
tkh (b) . Data obtained by the algorithm LPB for bottle-brush 
polymers with Nh — 515, and Nh — 35 are shown by dashed 
and dotted curves, respectively. 



FIG. 7: Log- log plot of reseated mean square height 
Rl{N, (t)/N^'' versus N (a) and ry = crN'' (b) with v = 0.588. 
Results are obtained for three choices of Nh and several 
choices of the grafting density a as indicated. Those unphys- 
ical data [Rh > 0.5Nb) due to the artifact of using periodic 
boundary condition are removed. The slope of the straight 
line corresponds to the scaling estimate from Eq. (|6|. 



write down the scaling ansatz in the thermodynamic limit 

as TV ^ oo [nj, 



Ri{N,a) = N'^R'{rj), rj = aN^ 



with 
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For checking this cross-over scaling ansatz, we plot the 
same data as shown in Fig. [TJ^a), but rescale the x-axis 
from N to T]. We see the nice data collapse. As r] in- 
creases, a cross-over from a 3D SAWs to a stretched side 
chain regime is indeed seen, but only rather weak stretch- 
ing of side chains is realized, which is different from the 
scaling prediction {Eq. (|6|}. In this log-log plot, the 
straight line gives the asymptotic behavior of the scaling 
prediction for very large r]. However, this is the first time 



we can see the cross-over behavior by computer simula- 
tions. This cross-over regime is far from reachable by 
experiments. On the other hand, it requires a lot of ef- 
fort to reach the regime where the theoretical prediction 
would apply, either the grafting density a has to be much 
higher or the side chain length N has to be much longer. 
The same situation is also observed as we check the 
scaling prediction for the radial distribution function [11] 
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Results are shown in Fig. [8] for bottle-brush polymers 
with side chain length N = 500 and A^ = 1500. As we 
keep all the grafting densities a fixed but increase the side 
chain length N^ one can see only in a rather tiny regime, 
the data seem to follow the same slope as predicted by 
the theory. Therefore, one might expect that finally the 
radial distribution function would follow the predicted 
scaling law as both a and N are very large. 
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FIG. 8: Radial distribution function p{r) plotted against r, 
for N = 500 (a) and iV = 1500 (b). 



For the second part of the simulations, let's first look 
at the snapshots of bottle-brush polymers under a good 
solvent condition, which contain N^ = 515 backbone 
monomers, A^ = (linear polymer), 6, 12, and 24 
monomers on each side chain, and the grafting density 
a = 1. As N increases, one sees that the correspond- 
ing conformations of bottle-brush polymers are rather 
different. The backbone becomes stiffer as the side chain 
length increases. This local intrinsic stiffness of the back- 
bone is quantitatively described by the persistence length 
£p. According to the scaling law of the mean square end- 
to-end distance of the backbone 17,8, 
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where i^ ~ 2.7 is the average bond lengths for the bond 
fluctuation model. Results of the rescaled mean square 
end-to-end distance of the backbone for TV = 0, 6, 12, 18, 
and 24 and for various numbers of backbone monomers 
Nij are shown in Fig. 10 ^a). The persistence length ip for 



fixed side chain length N is determined by the plateau 
for large A^5 since finally those curves of {Rl^)/{2£}jN^^) 
all show a smooth cross-over from a rod-like chain to a 
3D SAW. ip increases as N increases. 

The most common quantity to describe the structure of 



macromolecules is the structure factor S{q) for the whole 
bottle-brush polymer, which is estimated by taking the 
average of all independent configurations obtained from 
MC simulations, i.e. 



S{q) = 
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where c{fi) = 1 if f^ is occupied, otherwise c{fi) = 0. 
In experiments, S{q) can be measured by using static 
light scattering, small angle neutron scattering, and x- 
ray scattering, e.g. Ref. [19 . By choosing the accessible 
size of bottle-brush polymers for experiments, we have 
found the connections between our MC simulation results 
and the experimental data [7J. One example is shown in 
Fig. 10 ^b). As we normalize the structure factor S{q) -^ 



1 as g ^ 0, and rescale the wave factor q to qRg where 
Rg is the radius of gyration of the whole bottle-brush 

polymer, we see that the backbone length N^^^^^ = 400 in 
the experiment corresponds to the backbone length TV^ = 
259 in the simulation, and side chain length A/'^^^p^ = 62 
in the experiment corresponds to the side chain length 
A^ = 48 in the simulation. The grafting density a ~ 1 
for both cases. Immediately, we can translate that Inm 
~ 3.79 lattice spacings. 



IV. CONCLUSIONS 

In this paper, we study the bottle-brush polymers un- 
der good solvent conditions by using two kinds of lat- 
tice models, a simple coarse-grained model on the sim- 
ple cubic lattice and the bond fluctuation model. Due 
to the complex characteristics of bottle-brush polymers, 
we have proposed two algorithms, a variant of PERM 
and LPB depending on the interesting regime of length 
scales. With our extensive MC simulations, we show that 
the stretching of side chains in the interior of the bottle- 
brush polymer is weaker than the theoretical prediction. 
A convincing estimate of the persistence length ip which 
describes the intrinsic stiffness of bottle-brush polymers 
depending on the side chain length is given. We also 
give a direct comparison of the structure factors between 
our simulation results and the experimental data. The 
newly developed algorithm LPB has also been employed 
successfully to study the conformational change of bottle- 
brush polymers as they are adsorbed on a fiat solid sur- 
face by varying the attractive interaction between the 
monomers and the surface [20 . 
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FIG. 9: Snapshots of bottle-brush polymers with Nh — 515 backbone monomers, N side chain monomers, and with the grafting 
density a — 1. As N increases from (linear polymer chain) to 24, the backbone becomes stiffer. Quantitatively, the local 
intrinsic stiffness of the backbone is described by the persistence length ip as shown in Fig. 10 'a). 
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FIG. 10: (a) Reseated mean square end-to-end distance of the bottle-brush polymers, {Rlij)/{2ii)N^^), plotted against the 
number of monomers on the backbone, Nb, for the grafting density a — 1, and several choices of side chain length N. The 
corresponding persistence length ip for fixed side chain lengths N are given by the horizontal curves, (b) Structure factors S(q) 
plotted against the wave factor q. Simulation results are obtained for the bottle-brush polymer with Nb — 259, A/" = 48, and 
cr = 1 by LPB. Experimental data for the sample B2 with N^^""^^ = 400 and AT^^^p) = 62 are quoted from Ref. [19 . A snapshot 
of the bottle-brush polymer generated by LPB is also shown in (b) . 
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